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We represent the two - dimensional planar classical continuous Heisenberg spin 
P^ ■ model as a constrained Chern-Simons gauged nonlinear Schrodinger system. The 

hamiltonian structure of the model is studied, allowing the quantization of the 
theory by the gauge invariant approach. A preliminary study of the quantum states 
is displayed and several physical consequences in terms of anyons are discussed. 

1. Chern-Simons (CS) gauged models have recently attracted a great atten- 

tion in the description of several phenomena in condensed matter physics as the 
Fractional Quantum Hall Effect (FQHE) and the High Temperature Supercon- 
ductivity [1, 2]. In this context, models of the Landau- Ginzburg type have been 
proposed as effective field-theories [3, 4]. In the stationary limit, soliton solutions 
to the self-dual reductions of these models appear. These can be exploited to rep- 
resent Laughlin's quasiparticles and anyons [5]. On the other hand, the present 
authors showed that the self-dual CS system can be obtained, in the tangent 
space representation, from the stationary two-dimensional Landau - Lifshitz equa- 
tion (LLE) [6] . This is the classical continuous isotropic version of the Heisenberg 
model, obtained by the discrete quantum version in the framework of the coher- 
ent state representation [7]. In particular, soliton solutions of the former system 
turn out to be related to magnetic vortices of the LLE. More generally, the LLE 
in 2+1 dimensions takes the form of a constrained nonlinear Schrodinger system 
(cNLS) for two charged matter fields, coupled to a CS gauge potential [8]. This 
correspondence suggests an explanation to the phenemenological analogy between 
the behaviour of magnetic vortices [9] and Hall particles in an external magnetic 
field [10]. 

Sec. 2 is devoted to a brief review of the Landau - Lifshitz model in 2+1 



£ 



dimensions and its tangent space formulation, that is the CS gauged cNLS system. 
Furthermore, the conservation of some associated physical quantities is discussed. 
In Sec. 3 we study the lagrangian formulation and the hamiltonian structure of 
the CS gauged cNLS system. Finally, in Sec. 4 we apply the method of gauge 
invariant quantization to the above mentioned system. 

2. In 2+1 dimensions the LLE for the local magnetization S reads 

S t = SxV 2 S, (1) 

where S belongs to the 2-dimensional sphere S 2 . In terms of the stereographic 
projection 

S+ = S 1 +iS 2 = 2 Q -— 1 , S 3 = 1 ~ lC \ 2 (2) 

i + ICI i + ICI 

Eq. (1) becomes 

ji(t + (zz=2 — ^TT&G (z = x 1 +ix 2 , z = x t -ix 2 ). (3) 

With ferromagnetic boundary conditions, that is S(r,t) — ► (0,0, —1) for |r| — ► 00, 
the (xi,x 2 ) plane is compactified. Thus S defines maps of S 2 — + S 2 , which are 
partitioned into homotopy classes, labeled by the integer-valued topological charge 

1 f a ra c ., a c\ j2„ _ 1 f \Cz I - |C; 



Q = — S- (dtS x d 2 S) d l x = - / ^^ ^- d z x. (4) 
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The energy of the system takes the form 



1 f f^7C \ 2 j2„ _ A f\Cz\ \-\(z\ ,1, 



W = - / [VS ) d 2 x = A / lw ! ^ z ' d 2 x, (5) 

2j[ 1 J (l + IClf 

which compared with Eq. (4) leads to the Bogomol'nyi inequality W > 4n \ Q \. 
When the equality holds, the field £ has to be an (anti)holomorphic function of z, 
leading to the Belavin-Polyakov instantons [11]. 

In our approach the spin S is represented by the matrix S = S^a^ (a = 
1,2,3), where o~ a are the Pauli matrices. Then, we diagonalize S by resorting to 
a right-invariant local U(l) transformation g, such that S = gu^g -1 . The matrix 
g belongs to U(2) and has the form 

g = ±[2(l + S 3 )r 1/2 ( l +/ 3 ~f + S3 ) exp («r 3 |) , (6) 



where rj is an arbitrary real function. The corresponding chiral current J ^ (p = 
0, 1, 2) is decomposable into a diagonal and an off-diagonal part: 

/ M =,-^ = (i)^ + ( g ° 7»). ( 7) 

By definition, the chiral current J^ satifies the zero curvature condition 

d„J v - d v J^ + [J M , J v \ = 0. (8) 

Its off-diagonal and diagonal parts are 

D^q u = D v q^ (9) 

[D M , £>„] = - l - (5 M K - a^) = -2 (^ - ^„) , (10) 

where we have introduced the covariant derivative _D M = <9 M — f V^. 
In the variables g M , the equation of motion (1) becomes 

Qo = iD m q m , (11) 

which allows us to eliminate the function qo from Eqs. (10-11). Now, let us 
introduce the new fields 

iP±= {qi ^l q2 \ A = y -8(|^+| 2 + |^_| 2 ), A l = V l (i=l,2) (12) 
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and, correspondingly, the new covariant derivatives D^ = d^ — (i/2) A^. Then we 
notice that the expressions of ip± and Aj in terms of £ are 

^"2^ (i + ICI 2 ) ' J (l + ICI 2 ) J ^' { } 

For |U = 1 and v = 2, we see that Eq. (9) takes the form 

1 = {D 1 + iD 2 ) ^_ - (D 1 - iD 2 ) ^+ = 0. (14) 

This equation can be interpreted as a geometrical constraint for the CS gauged 
NLS system, arising from the remaining equations (9) and Eq. (10), i.e. 

iD ^± + (D t 2 + D 2 2 ) ij>± + 8tt |V>± | 2 V± = 0, (15.a) 

d 2 A 1 - d x A 2 = 8tt ( J+ - J ") , (15.6) 

d Q Aj - djA = -8tt £ji ( J+ - J,~) , (15.c) 

where 

4 = i^±i 2 5 j ; ± =i(^±A^±-^±A^±) a = 1,2). (i6) 



Here we introduce the charge density currents satisfying the continuity equation 
<9 p J^ = 0. Then, the matter fields ip± carry a conserved electric charge and we 
can interprete the quantity 



Q = J (J + - Jo) J 



X 



(17) 



as the (conserved) total electric charge. In terms of tp± it is easy to see that Q 
and Q coincide, thus we are led to the quantization of the total electric charge. 
Furthermore, Eq. (15. b) implies the quantization of the magnetic flux: $ = 



jBd 2 x = J(d 1 A 2 -d 2 A 1 )d : 



x 



-8tvQ. 



Moreover, the continuity equations for the currents J^ predict that the total 
number of particles (anyons) 



Af 



J ( 4 + j -) d< 



X 



(18) 



is also a conserved quantity. Now, comparing Eq. (5) with (18) , we obtain W = 
AtiM. Then, in the tangent space formulation, the Bogomol'nyi inequality takes 
the form \Q | < A/". Of course, the equality holds when tp + , or ift_, vanishes. In 
such a case, the system (14 - 15) reduces to the self-dual CS model [2]. This model 
is completely determined in terms of the solutions of the Liouville equation. This 
simple observation enables us to identify the magnetic vortices, or equivalently the 
instantons of Belavin-Polyakov type, with the solitons of the self-dual CS model, 
that is with the anyons. In particular, from (17) we have an explicit correspondence 
between the magnetic topological charge and the anyon electric charge, whose 
quantization finds here a geometrical interpretation. In ref. [8] the more general 
static case (Di + iD 2 ) i/i- = (-Di — iD 2 ) V>+ = is treated. This leads to a 
special reduction of the euclidean conformal affine Toda system, which contains 
the sinh-Gordon, or the Poisson-Boltzmann equation, in a particular gauge. 

Finally, we notice that the constraint (14) is identically satisfied using the 
expressions (13). Indeed, we find the identity 



(D 1 + iD 2 ) ^_ = (D 1 - iD 2 ) ip+ = \ 



S 2 2 
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Equating to zero the quantity in the square brackets, we recognize the stationary 
LLE in the stereographic variable (see Eq. (3)). Thus, all static solutions to Eq. 
(1) are mapped, modulo gauge transformations, into solutions to the sinh-Gordon 
reduction studied in ref. [8]. 

3. System (15) can be derived from the Lagrangian 



C= l - {^D^ - A)« - D a ^D a ^ + 4n (^a 3 ^) (^) + -^- e^ A A^d v A x , (20) 
A 6Ztv 



where we have used the spinor notation ifj T = (ip+,ip_), with ip = ifj + o~3, and the 
summation over a = 1, 2 is understood. Let us observe that the strength of the self- 
interaction and the CS coupling are fixed by the original spin model. The specific 
ratio between these constants allows us to perform the self-dual reduction, which 
is a very special subcase in more general and the phenomenological theories of CS 
type [2]. Moreover, in the effective theories for FQHE the CS coupling is related 
to the so-called filling factor v, which in the present case is 1/2, where the units 
fi — c — 1 and e = 1/2 have been used. This result is consistent with the so-called 
chiral spin fluid model [1 - 12]. Moreover, recent experimental observations in 
double-layer samples indicate the compatibility of u, with values which are inverse 
of even numbers [13]. 

A final remark concerns the parity invariance of the Lagrangian (20), which 
changes the total sign under the transformation 

^ + ^V-, D 1 + iD 2 <-> D 1 - iD 2 . ^ ' 

So that equations of motion and the constraint (14) are both left invariant. This 
fact is an obvious consequence of the parity invariance of the spin model (1) and 
it is another example, against the common belief, that the presence of the CS 
interaction breaks the P and T invariance. 

Now let us enter into the detail of the hamiltonian structure of the model 
(14 - 15) by resorting to the first order Lagrangian method proposed in ref. [14], 
which implements the classical Dirac method for constrained Hamiltonian systems 
[15]. 

First, we consider as dynamical variables £ = [ip+, ip+, ip~, ip-,Ai,A2) and, 
defining a - (§ $+, -\ ip+, -§ V>_, \ ip_, ^ A 2 , -g|^ Ai), we rearrange (20) in 
the form 

£' = a m e m -ft(£) + ^4 ri, (22) 

H = D a ijD a ij - An (^a^) (« (23) 

Fl = JIT = W + F" e * id * A i = °- ( 24 ) 

Thus CJ describes the system (15 a,c) constrained by Eq. (24), that is the Chern- 
Simons version of the Gauss law (15. b). In (22) Aq plays the role of Lagrange 
multiplier. Furthermore, 7i would be the density of the Hamiltonian H . In fact, 
the corresponding Euler-Lagrange equations take the hamiltonian form 

e = fi J - 1 ^ = {e,H}, (25) 

where f m i = dai/d^ m — da m /d£ is a nonsingular symplectic matrix and we have 
defined the Poisson brackets by {£ m , £*} = /~ ; S 2 (x — y). Specifically, in our case 
the nonvanishing Poisson brackets are 

{V>±(x),^±(y)} = Ti 5 2 (x-y), {^(x),^(y)} = lQire^S 2 ^ - y). (26) 



By using Eq. (33), one easily shows that f i = {Ti, H} = 0; consequently, in phase 
space the Gauss' law defines an invariant submanifold under the dynamics. On 
the other hand, Eq. (14) provides the further complex primary constraint 7 = 0, 
whose real and immaginary parts will be denoted by 1?2 and 1?3. By evaluating 
the Poisson brackets, we find that {7, H} ps 0, i.e. 7 is vanishing only on the 
submanifold defined by (14) or, in other words, Ti and 1?2 are weakly invariant 
under time evolution. Furthermore, they do not produce secondary constraints. 
However, since {Ti, 1?2, Ts} is a set of first class constraints (i.e., {Tj, H} ps and 
det({rj,rj}) ps ), we need a gauge-fixing condition, for example the Coulomb 
gauge condition 1?4 = d{A l ps 0, which is the most used in the litterature [16, 
17, 18]. So, we are able to compute the Dirac brackets, defined by {£ m ,£'}.D = 
{f",^} - {r\r a }(0 M ,{r a ,,£'}, where d ss , = {r s ,IV}. In our case, the non- 
vanishing matrix elements d ss i are distributions in terms of 5 2 (x — y) and its 
derivatives. Then, inversion can be performed only in implicit form. In order to 
avoid this type of problem, one can look for a different gauge-fixing condition. For 
instance, a more general non-local gauge fixing can be choosen [19], or the so-called 
" axial gauge" [20] . Nevertheless, in all generality the corresponding Dirac brackets 
are rational in the field variables. Then a quantization of the model seems to be 
rather difficult, because of the problem arising in the ordering and the inversion of 
the operators involved. This is one of the main reasons which address us to look 
for an alternative approach, which is invariant under the choice of the gauge. 
4. In this Section, we will adopt a slightly different point of view, in order to 

tackle the above mentioned problems. Moreover, we would like to understand the 
quantum meaning of the geometrical constraints 7, as seen separately from the 
basic CS-NLS field theory. Then, a possible route consists in quantizing system 
(14 - 15) in two steps: 

1. first, we apply the procedure of the gauge-invariant quantization to system (15) 
following the approach suggested in ref. [21]; 

2. then, we use the constraint 7 to select the physical states among all the gauge 
invariant quantum states. 

First of all, let us change the hamiltonian description by taking 

w' = w(0 + /or + ^4ori, (27) 

where we have as primary constraint To = tcq = 0, ttq being the momentum 
conjugated to the canonical variable Aq, such that {Aq,ttq} = 5 2 (x — y). This 
constraint is implicit in the previous formulation, since dC/dAo = 0. Ti is now 
a secondary constraint, arising from the Poisson bracket {no,H}. Finally, in the 
expression (27) /o is an arbitrary function, which defines the evolution of Aq. 
Furthermore, noticing that the special canonical pair (A\, A2) breaks the rotational 
covariance of the theory, it is more convenient to consider the decomposition of 
the vector potential into a longitudinal and a transverse part Ai(x) = A\ + Aj = 
2diT]('x) + eij9jx(x). Recalling that the magnetic field is B = e^diAj = — <9 2 x, we 
derive formally the useful relation 

^(x) = -e, J (a- 1 S)(x), (28) 
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where we have defined d- 1 / (x) = ^rd Jin |x — y\f(y)d 2 y. Analogously, it is 
convenient to express ip± in terms of the canonical variables (P±, Q±) by 

il>± = ^={P±TiQ±). (29) 

Hence, the Poisson structure associated with the hamiltonian system given by (27) 

is 

(Q±(x), P±(y)} = 5 2 (x - y), {A (x), 7r (y)} = 5 2 (x - y), 

{ J B(y),r 7 (x)} = 87r5 2 (x-y). ^ 

Furthermore, the equations of motion are supplemented by the first class con- 
straints ITo ~ Ti pa 0. 

Now, we are look for a suitable canonical transformation, such that some 
of the new momenta are equal to the constraints. In such a way, the coordinates 
canonically conjugated to these momenta have an arbitrary time evolution and are 
remnants of the gauge invariance of the theory. The remaining canonical variables 
define the prominent dynamical degrees of freedom of the theory. 

First, let us denote by (Q±, P±), (Aq, tTq), (rj*, ir^) new canonically conjugated 
coordinates, where the momenta it* have to be set equal to the constraints I\. This 
can be done through a suitable generating function F = F (Q±, Aq, 77; P±, 7Tq, nf) , 
satisfying the set of generalized Hamilton- Jacobi equations 

* r (n A 9F dF dF\ 

The integrability of this system is assured by {Fo,Fi} = 0, moreover its general 
solution is given by 



F = n A + ^ + -P + 2 - -Pjj V + 

J ^P; 2 - Q\dQ + + J yJp?-Q*dQ-. 

The corresponding canonical transformation reads 

tto = n* , A = A* , -±B = n$ + \Pf - \P1\ 77 = 77* 
Q± = PI sin ( §| T V*) , P± = P± cos ( §1 T V* 

Finally, the new gauge-invariant degrees of freedom are 



(32) 



(33) 



1 / Q* 

$ ± = — P|exp( T i-pf ) =^±exp(-i?7), ir*± = $±, (34) 

which fulfill the canonical brackets 



{$ ± (x),$ ± (y)} = T z5 2 (x-y). (35) 

In terms of these new variables, the Hamiltonian density becomes 

n = (d a + iAT)* + {d a - |^>+ " (d a + ^)$_(0 a - iAl)<S>_ 

(36) 

-4 7r(|$ + | 4 -|$_| 4 )]+/°7r* + iA*7r*, 

where A^ (x) = 8iv e a bd^ (l^+l 2 — l^-| 2 )( x )- Furthermore, the constraint 7 takes 
the form 



-, = (d + --Al)S:-(d---A T _)<t + , (37) 



where d± = di ± i 82 and A± = Af ± i A%. 

Now, let us observe that the relations (13) are algebraically invertible in terms 
of gauge invariant quantities when the topological density is non vanishing, namely 

C= 4^ |$ + | 2 -|$_| 2 • (38) 

Then, we can directly reconstruct the spin field in terms of $±, although it turns 
out to be a non local field in view of the expressions for A±. In this regard, let 
us observe that, in the one-dimensional completely integrable case, the relation 
between £ and the NLS matter field is also analytic in nature, but we need to 
solve the associated linear spectral problem. 

Furthermore, keeping in mind that the current densities (16) are manifestly 
gauge invariant, we can give the expressions of the linear and angular momentum 
for the model (14 - 15) in a gauge invariant form. 

Firstly, in opposition to what happens for the electric charge densities, the 
linear momenta P^ = J J^ d 2 x, associated with the <fr± fields separately, are no 
longer conserved, but they evolve according to Newton's law 



dt 



"m ~~ ° e mn I ("A) ^n ^0 ^ n I" 1 X - \^) 



We notice that both <&± are subjected to the same force. Therefore, the total 
kinetic momentum 

p m = p+-p- = J ( $+a m $+ - $_a m $_) d 2 x + c.c. (40) 

is conserved, the " — " comes from the "negative mass" associated to $_ . The 
total angular momentum M = M + — M_ is also conserved , with 



M±= (xxJf - x 2 jf) d 2 x = £± + s± 



(41) 



and 

l± — —i f$>± (xid 2 - x 2 di) §±d 2 x + c.c, , ■. 

s± = f\$±\ 2 (AT x 2 - A^x 1 )d 2 x. ^ ' 

The quantities I — 1+ — I- and s = s + — S- are interpreted as orbital and spin an- 
gular momentum, respectively. In particular, the general result s = 2 Q 2 suggests 
that vortices have non zero momentum also in the static case. 

Finally, since the magnetization M. = J S d 2 x is a conserved non-local quan- 
tity for the LLE system, via Eqs. (38) and (2) we find a further gauge invariant 
conserved quantity for the system (14 - 15) , related to the original global 577(2) 
symmetry. 

We are now ready to pass on the quantum theory in a straightforward manner, 
just by replacing the canonical brackets by i times the commutators (or the anti- 
commutators, if the matter fields are supposed to be fermionic) . At the quantum 
level, the first class constraints r$ become the operators I\, which must annihilate 
the physical states. This fact implies that the quantum states are independent of 
Aq and invariant under time-independent gauge-transformations. Therefore, the 
gauge - invariant operators of the theory must commute with IV Then, recalling 
the Poisson brackets for the Hamiltonian system described by (36), we can estab- 
lish the correspondence <fr± — ► <5± and $± — ► $j_, with the equal time (anti-) 
commutation relations 

[i>±(x),i4(y)] ± = ±5 2 (x-y). (43) 

In the physical gauge-invariant subspace of the full Hilbert space, the Hamiltonian 
is given by 



H = J{-& + {d a ~ ^) 2 ^ + + *!(0 O - | Al) 2 ®. 
-4tt : ($^$+) 2 - ($!i>_) 2 :}d 2 x, 



(44) 



where the normal product operator : : and A^ = 8ne a bd b 1 ($] ) _l ) + — $Ll ) _)(x) 

have been used. The fields $± are local and invariant under local gauge transfor- 
mations generated by the first-class constraints. Therefore, they may be exploited 
to study the breaking of global simmetries. 

The quantum version of relations (30) leads to [i?(x),r)(y)] = 87ri5 2 (x — y), 
which implies the relation 

exp (i fj(y))B (x) exp (-i fj (y)) = B (x) + 8 n 5 2 (x - y) . (45) 

This result is used to prove that <&± are the gauge invariant operators which create 
a charge-solenoid composite, having magnetic flux equal to ±8 7r . 

The Hamiltonian (44) commutes with the particle number operators 

N±= [ & ± $ ± d 2 x. (46) 



Hence, we can diagonalize them simultaneously via the Fock representation for a 
quantum state with N = N + + iV_ particles, which is given by 



\N+, N_ > = / Jl d 2 x+ Y[ d 2 Xj * (x+, . . . , x+ + , Xl , . . . , x" 

J i=l j = l 

*t (x+) ...*t_ ( X +J <|t_ ( xr ) ...6L (x^) |0 > 



(47) 



where the vacuum state is defined by $±|0 >= 0. The function \1/ is an arbitrary 
element of the Hilbert space L2{1Z 2 ( N++N ~>] obeying the Schrodinger equation for 
(iV+ + iV_)-bodies. Furthermore, the state (47) is gauge invariant by construction. 
Thus, we have selected in the full Hilbert space only those states which contain 
the whole physical information. 

Nevertheless, up to now we have dealt with a quite general "anyonic" field 
theory, without using the constraint 7, which in some sense takes into account the 
properties of the original spin field model. As a consequence, "quantizing" Eq. 
(37) in the form 

7 = £>+$_ (x) -£>_$+ (x) , (48) 

we should obtain the physical states by requiring that / y\N + , N_ >= 0. However, 
for arbitrary occupation numbers (N + ,N_) this equation leads to complicated 
relations for the wave function \1/. Here we display the result only for the finite 
boson case |iV + ,0), namely 

vi>(x+,---,x+ + )=^^ 1 +,---,^+_)nl x ^- x J + ~ 2 ' ( 49 ) 



+ , 

i<j 



where JF is a holomorphic function of its arguments. By using a singular gauge 
transformation [23-24], the expression (49) takes the form of the Laughlin wave 
function [25] , describing a condensate state of bosonic solitons, which may be re- 
lated with quantum disorderd state of the original ferromagnets [12]. Furthermore, 
it should be proved that the quantum state corresponding to the wave function 
(49) is actually an eigenvector of the Hamiltonian (44) . This work is in progress, 
and we have strong indications in the positive sense, because at the classical level 
the Hamiltonian (23) can be written as a linear combination of the constraint 7 
and Ti. 
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